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ABSTRACT 

We test the applicability of the Gallavotti-Cohen fluctuation formula on a nonequilibrium version of the periodic 
Ehrenfest wind-tree model. This is a one-particle system whose dynamics is rather complex (e.g. it appears to be 
diffusive at equilibrium), but its Lyapunov exponents are nonpositive. For small applied field, the system exhibits 
a very long transient, during which the dynamics is roughly chaotic, followed by asymptotic collapse on a periodic 
orbit. During the transient, the dynamics is diffusive, and the fluctuations of the current are found to be in agreement 
with the fluctuation formula, despite the lack of real hyperbolicity. These results also constitute an example which 
manifests the difference between the fluctuation formula and the Evans-Searles identity. 

I. INTRODUCTION 

In molecular dynamics simulations of fluids in nonequilibrium stationary states, Evans, Cohen and Morriss 
1^] discovered a remarkable relation for the fluctuations of the entropy production rate. This relation links 
in a striking fashion the microscopic reversible dynamics of certain particle systems in a stationary state, 
to the corresponding irreversible macroscopic dynamics. Inspired by these flndings, Gallavotti and Cohen 
proved ||^ the fluctuation relation for a wide class of systems, directly from the dynamics of their constituent 
particles. Their proof was based on the following: 

Chaotic Hypothesis (CH): A reversible many-particle system in a stationary state can be regarded as a 
transitive Anosov system for the purpose of computing its macroscopic properties. 

The ensuing result is now known as the Gallavotti-Cohen Fluctuation Theorem (GCFT). Near equilibrium, 
the GCFT implies both the Onsager and Einstein relations |^ and can therefore be interpreted as an 
extension of them to far-from-equilibrium situations. 
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On the other hand, quoting from the review article ||^: "... in concrete cases not only it is not known 
whether the system is Anosov but, in fact, it is usually clear that it is not ... Hence the test is necessary to 
check the CH which says that the failure of the Anosov property should be irrelevant for practical purposes. " 
For "practical purposes" means that the calculation of quantities of physical interest is not affected by the 
deviations of the dynamics from the ideal case of an Anosov flow. Therefore, numerical or real experiments 
are required to test the applicability of the CH, and to identify its range of validity. Several papers have 
been devoted to this purpose (see, e.g. Refs. (5|-||]), while other papers have investigated the possibility of 
observing fluctuation relations similar to that of the GCFT in different contexts (see e.g. |ic|-[l^]). These 
works show that the CH is appropriate in the interpretation of: (a) the numerical results obtained for 
two-dimensional systems of hard-core particles [^|j8j; (b) some experiment on liquids undergoing Benard 
convection |^ ; (c) numerical simulations of two-dimensional turbulent fluids |^ ; (d) the heat transport along 
chains of anharmonic oscillators |^]. At the same time, the papers |l|Jl| extend the validity of the GCFT 
to stochastic dynamics, including rather general Markov processes. 

^From all the mentioned examples, one can indeed conclude that -according to the original intuition of 
Ref. - the CH effectively works for a definitely wider class of systems than that of topologically mixing 
Anosov diffeomorphisms or flows. Of course, despite of the lack of uniform hyperbolicity (or smoothness, or 
both) all these time-reversal invariant systems share the common property of being strongly chaotic, in the 
sense that they have (possibly many) positive Lyapunov exponents. ^ 

One question comes to the fore: given a time-reversible dynamical system, what is the minimal degree 
of "complexity" required for its microscopic dynamics to verify the fluctuation relation? Alternatively, one 
may ask how "weakly chaotic" can be a system which verifies the CH. 

These are rather natural questions in the framework of statistical mechanics, where similar problems 
have traditionally been considered. For instance, the assumptions of ergodicity or of molecular chaos, are 
universally accepted for equilibrium systems, in spite of the well known fact that such properties are not 
only exceedingly difficult to prove in practical cases, but are violated in most of them. Nevertheless, these 
assumptions lead to the correct physical predictions, and provide a mechanical foundation to thermodynamics 
by linking the latter to the microscopic dynamics. The reasons of their success are hidden in the interplay of 
extremely different time and length scales, and in the large numbers of particles which constitute macroscopic 
systems (see, e.g., P,|T^-p^ for a discussion of these topics, and Ref. for a recent work on the role of 
different time scales in classical gases). 

In the present paper we approach those issues by studying a modified version of the Ehrenfest wind-tree 
model. It consists of a particle bouncing elastically in an array of fixed polygonal scatterers. The original 
version (randomly distributed and square obstacles) was introduced to study the validity of Boltzmann 
equation and has been very recently reconsidered |p^ . The most remarkable properties for our purposes 
is the fact that flat boundaries prevent chaotic behaviour as no defocussing of nearby trajectories occours. 
Nonetheless, the model has "good" statistical properties: the iJ-thcorem holds, and the moving particle 
gives rise to a true brownian motion. 

What will be considered here is a nonequilibrium version of the model, with an external field and a 
Gaussian thermostat (see, for the details, the next section). The behavior of the system for nonvanishing 
field is essentially the following: asymptotically the dynamics is trivial, namely the particle collapses onto a 
periodic orbit. Before this, however, there is a long transient (of duration growing as e~^, if e denotes the 
intensity of the applied field) during which the motion is basically diffusive as in the equilibrium case. Our 
main result is that the GCFT holds during such a transient. One may therefore speculate that, although 
chaoticity is in principle necessary to guarantee the validity of the GCFT for all times, the latter may still 



'Strictly speaking, such remark applies only to deterministic dynamics: nonetheless, stochastic systems are consistent 
with the presence of chaotic underlying dynamics 



2 



retain its meaning even in the absence of chaoticity, for trajectories of large ahhough finite length. 

Beyond this, our results also contribute, in our opinion, to clarify questions which have been recently 
raised in the literature. First of all, they provide an example in which the different origin of the GCFT and 
of a previous result, known as the Evans-Searles Identity are evident (cf. Section III). Moreover, our 
findings give further support to the claim that observing diffusive behaviour does not constitute by itself 
a proof of the chaoticity of the microscopic dynamics (cf. Section II). 



II. THE MODIFIED EHRENFEST GAS 



The periodic Ehrenfest gas consists of one moving particle of mass m, which is elastically scattered by a 
set of rhomboidal fixed obstacles arranged on a triangular lattice (see fig. |l|). Between any two collisions, the 
particle moves freely along straight lines. The periodic structure of the lattice allows us to follow the motion 
of the particle, looking at the periodic image of its trajectory in just one hexagonal cell: the elementary cell 
evidenced in fig. ||. As the sides of the scatterers are flat, the dynamics cannot be chaotic: all of the the 
four Lyapunov exponents vanish. Nevertheless, for generic (i.e. irrational) values of the internal angles of 
the rhombus, the system is expected to be ergodic, so that in particular a long enough trajectory fills the 
available phase space (see for a review on the subject, and for more recent results). 

Similarly to the case of the nonequilibrium Lorentz gas p^, we modify this model by adding an external 
field of intensity e, and introduce a Gaussian thermostat, which constraints the kinetic energy of the particle 
to its initial value K. Let (x, y) be the position of the particle, and (j)x,Py) be its momentum. We fix m = 1, 
K = 1/2 and let the field point in the positive x-direction. The equations of motion for the free flights thus 
read 

'x^Px\ px = -apx+e 

with a — epx ■ (1) 

y=Py] Py^ -apy 

The effect of the external force is that the segments of trajectory between subsequent collisions get curved, 
but because of the simple form of Eqs.(|l|) they can be computed analytically (see the Appendix). 

We performed numerical simulations of the model for small and moderate fields, 10^"' < £ < 1 by evolving 
eqs. (If) starting from generic initial conditions. The asymptotic state was always found to be a periodic 
orbit. J Accordingly, upon switching on the field, two of the four Lyapunov exponents remain zero (the 
exponents corresponding respectively to the conserved kinetic energy, Ai, and to the direction of the flow, 
A2) while the other two are numerically seen to approach a negative or vanishing value: 

= Ai = A2 > A3 > A4 . (2) 

In practice, for small fields like those mostly considered here, the exponent A3 is often so small that numer- 
ically we can hardly distinguish it from zero, while A4 takes definitely negative values. 

The mere existence of a trivial asymptotic motion for e =/= does not however exclude quite complicated 
behaviour. Indeed, as dissipation is weak for small e, a long transient is required to reach the periodic orbit. 
During the transient the motion of the particle looks rather erratic, and covers almost uniformly a large 
fraction of the phase space fl. Furthermore, the behaviour of the system on this time scale appears to be 
almost stationary, and can be described in a statistical way. 

To illustrate these facts, and visualize the dynamics, it is convenient to introduce the usual "bounce map" 
of billiards. Precisely, each collision is assigned coordinates (s, cos^p), where s is the distance of the collision 



^We recall that periodic orbits in the elementary cell may be open in full phase space, with a total displacement of 
an integer number of lattice vectors for each period. 
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point from, say, the rightmost angle of the scatterer, measured along its perimeter, and ip is the angle 
between the outcoming momentum and the side of the rhombus (positively oriented in the counterclockwise 
direction); the billiard dynamics then defines a map 



B : (s„,cos-0„) 



I — !• 



(s„+l,COS V'n+l) 



(3) 



where (s„,cosV'n) denote the coordinates of the n-th collision. For e — 0, the map turns out to be area 
preserving, so that ergodicity corresponds to uniform filling of the square (or, better, of the cylinder) [0, L) x 
[—1,1], L denoting the overall lenghth of the border. 

As remarked in the Introduction, the asymptotic behavior of the system is trivial, namely all trajectories 
eventually approach a periodic orbit. However, before reaching the asymptotic regime, the system exhibits 
a long transient, up to some number of collisions Uc (see later for an appropriate definition of ric) during 
which the dynamics looks nontrivial, as if the system were chaotic. The two regimes are illustrated in Fig. 
|[ for case e = 0.01, which has ric ~ 3 x lO'*. The left and right panel of the figure report, respectively, the 
first and the last 5,000 iterates of the map, out of a trajectory of 10^ collisions. Clearly, during the transient, 
5,000 iterates are sufficient to roughly cover the square, while asymptotically one is left with few isolated 
points. A closer inspection, by means of hystograms of the density of points in the square (not reported 
here) shows that, for smaller s and correspondingly larger Uc, the iterates of the map during the transient 
fill more and more uniformly the square. 

A more quantitative characterization of the two dynamical regimes is achieved by considering a large set 
of initial data {(a;o*\yo*'''Pio'^'yo)}i^i' picked up at random with uniform distribution in the phase space, 
and measuring the variance tr^ = {{xn — a^o)^ + iVn — yo)^)i where {xmUn) is the actual position of the 
particle at the n— th collision, and (.) denotes averaging over initial data. This corresponds to an ensemble 
average, for a non-interacting gas of indipcndentely thermostatted particles. As is clear from Fig. ^, one finds 
a crossover, at rt w iic, from diffusive (cr^j ^ n) to ballistic (cr^ n^) behaviour, the latter corresponding to 
the approaching asymptotic periodic orbit. This is in fact our definition of Uc- By varying e, one finds that 
ric is inversely proportional to e^, see Fig. ^. The divergence of the transient for e ^ is fully consistent 
with the conclusions of |l^ ] , according to which genuine diffusive behaviour is found for the wind-tree model 
without an applied field. As already mentioned in the Introduction, this seems to contradict, or at least to 
weaken, the claim (see pl[|) that observing diffusive behaviour in a given physical system is a good indicator 
that the corresponding molecular dynamics is chaotic. 

Similarly to the case of the Lorentz gas, the presence of the field induces also here an average drift of the 
particle along the field direction, so we can define the "current" in the system as 



i.e. as the ensemble average of the time average of the component of the particle momentum along the field 
direction. All the cases we considered displayed a well-defined positive current {px)t > 0, both for moderately 
long t (in the quasi-stationary transient state) and for very long t (in the asymptotic state). 



The discussion of the previous section can be summarized by saying that, for generic initial conditions, the 
particle spends a considerably long time exploring practically all the available phase space, until it eventually 
reaches the periodic orbit. Moreover, the time duration of the transient Tc — toHc, tq being the mean flight 
time, diverges as upon decreasing the field. So, if we restrict the attention to time scales shorter than 
Tc, it makes sense to study the statistical properties of the fluctuations of a given observable, and compare 
the result with the prediction of the GCFT. To this purpose, we first need to adapt the latter to the present 




(4) 



III. THE FLUCTUATION RELATION 
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case. Let T > 0, M E IN and take r — T/M. In our framework, the GCFT may be replaced by the following 
conjecture, based on our numerical observations: 

Conjecture: Consider a periodic billiard with flat scatterers, which is ergodic at equilibrium. Let the 
particles be subject to external driving and to a Gaussian thermostat. Then, there is a critical time Tc such 
that for T <^ Tc, and for sufficiently large M and t, the following holds: 

In / ' = ez + o{ez) , (5) 



where tx^ denotes the probability distribution of the quantity 

^ ^ {Px)t 

{Px)t 



(6) 



computed by subdividing a simulation of length T in segments of length t, and by recording the observed 
frequency of occurrence of the values z. 

In the original formulation of Ref. the fluctuation relation holds in the limit of large T and r, for fixed 
£. Instead equation makes sense only for T and t finite. Alternatively, the limits of large times and of 
small e should be taken simultaneously, with T well within the transient Tc ^ £~^, when the observables are 
actually fluctuating. The correction term o{ez) accounts for the observation that reducing e at flxed z, or 
reducing z at flxed s, the left hand side of eq. is well approximated by ez itself. 

Under the above limitations, we checked the validity of Eq.(^) for several values of e. In the numerical 
computations, it would be convenient to consider small fields as they correspond to larger Tc and entail better 
statistics for z on each run. On the other hand, the accuracy of the simulations worsens with decreasing field, 
as usual in nonequilbrium molecular dynamics (see, e.g., jl^), thus a compromise between these contrasting 
needs is required. The chosen e values range between 10"'^ and 10^^ and are given in the figures along with 
the values of t. Incidentally, note that these values suffice to go well beyond the linear regime of irreversible 
thermodynamics, which for models of this class can be estimated to correspond to fields of 10~^ or smaller 
(cf. Ref. ||l6|, point 2 of Discussion). 

As usually done for billiards we decided to cut our long trajectories in segments of a fixed number of 
collisions, rather than of a given duration. Hence, the values of r reported in the figures are to be understood 
as the average times necessary to undergo rir subsequent collisions with the scatterers. At variance with 
||], but similarly to we did not decorrelate the successive trajectory segments in order to have better 
statistics. A further average over an ensemble of iV ~ 10^ independent trajectories was also performed. 

Some typical results, reported in Fig. ^, seem to vindicate the above conjecture. In particular, the 
agreement with eq.(^ improves with increasing t, and the range of validity of the formula gets wider by 
diminishing the field. Fig. |^ illustrates the non-gaussian nature of the probability distribution for the 
same values of the parameters. Indeed, tTt-{z) appears to be slightly asymmetric and to have approximately 
exponential tails. 

The conclusion we draw is that, in spite of the absence of any source of chaoticity, the complex (diffusive) 
dynamics of our system in the transient states is almost indistinguishable from genuine chaotic dynamics. 
One may therefore wonder about the origin of the randomness in the wind-tree. A possible explanation is 
the following: it frequently happens that out of two nearby trajectories, one collides with a scatterer, while 
the other just passes by, without being scattered. This results in a separation of orbits for short times until, 
eventually, dissipation takes over and both orbits converge to the asymptotic one. This effect can be easily 
probed by computing the finite-time Lyapunov exponent A3(n), for n <^ Uc- The latter measures precisely 
the average separation of close enough initial points after a fixed number of collisions. The results reported 
in Table 1 show that the time scales over which the finite-time exponent appears to be positive covers all 
the time range of validity of the fluctuation relation (cf. Fig. ^). 
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n 


e = 0.01 


e = 0.005 


4 


3.62 10"^ 


3.60 10"^ 


9 


2.13 10"^ 


2.14 10"^ 


19 


1.19 10"^ 


1.22 10"^ 


39 


6.96 10"^ 


6.87 10"^ 


79 


4.26 10~^ 


3.96 10"^ 


159 


2.99 10"^ 


2.39 10"^ 


319 


2.44 10"^ 


1.67 10"^ 


639 


2.01 10"^ 


1.40 10"^ 



TABLE I. Finite time Lyapunov exponent \3{n), for different number of collisions n. Tfie ensemble size is A'^ = 10^. 
The trend for larger ensembles is similar. We observe that only for ric, ^^{n) will change sign. 



IV. CONCLUDING REMARKS 



Our analysis indicates that the consequences of the CH are valid (albeit in a restricted sense) for non- 
chaotic and reversible particle systems, provided that their trajectories are sufficiently "unpredictable" for 
relatively long times. In this respect, the difference between real chaotic systems and our wind-tree seems 
only to consist of the possibility to extend to infinite times the validity of the CH and of its consequences. 
Similar limitations on observation times are usual in statistical mechanics. 

Admittedly, the example discussed here is not completely generic at least for two reasons. First, it is 
fairly artificial to start with a system that is ergodic at equilibrium without being actually chaotic (even in 
a weak sense). Second, although other models similar to ours can be conceived, in general one expects the 
asymptotic state to be chaotic, even in the presence of external fields. Nevertheless, our model is important, 
in our opinion, from a theoretical point of view, because it provides a limiting case in which the consequences 
of the CH can still be applied. Accordingly, this also suggests that the CH can be successfully applied to 
systems with slow decay of correlations or long-time tails, strongly deviating from the true Anosov systems. 

As a final remark, we observe that our results are perhaps of some interest in the present debate on the 
relation between the GCFT, and one identity previously obtained by Evans and Searles (ESI) ||l^,^ . The 
ESI concerns any time reversible dynamical systems, like ours, and the Liouville mesaure /il on the phase 
space Q. of such systems. In particular, let Ep C he the subset of initial conditions of trajectories along 
which the phase space contraction is e"^^"^"^ after a time T > 0, where (.) represents a (stationary state) 
average. Then, the ESI can be expressed as pO[ |: 



^J'L{E-p) 



eP(-)^ . (7) 



This equation is formally similar to the fluctuation relation of the GCFT, especially if one takes very large 
T, so that the computed quantities characterize the stationary state of the system. However, the GCFT and 
the ESI cannot possibly refer to the same physical quantities. Indeed, if T is large in the modified Ehrenfest 
gas, there are no fluctuations of the phase space contraction at all. Accordingly, as observed above, the 
fluctuation relation does not apply. On the contrary, the ESI (0) retains its meaning, showing that the 
quantity on the left hand side of the equation cannot be interpreted as probability of flucutations. The 
ESI, instead, correctly gives a relation for the probabihties of "trajectory hystories" with opposite phase 
space contractions. 
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APPENDIX: LYAPUNOV EXPONENTS 

In this appendix we give some expressions for the (finite time) Lyapunov exponents of the thcrmostatted 
Ehrenfest gas. Let us consider the coordinates x, y, and 9, where the latter is the angle formed by the 
momentum vector {j)x,Py) with the x axis. Because the kinetic energy of our system is constant, K = 1/2, 
the coordinates {x,y,6) suffice to describe the dynamics which, going from one coUision to the next, can be 
split in two stages: 




.y'o, 




(8) 



Here F, the free flight between two successive obstacles, is explicitly given by (see eqs. (3) and (4) in p^ ) 

tan ^ = tan y g-'^^^*"'^"'^") (9) 

, 1 , sin6'n , , 

x'o^xo--\n 10 
e sm 6*0 

J/o = ?/o (11) 

e 

where r, the flight time, depends on the initial point (^Oi Vo)- In turn, the map C represents the collision 
with the scatterer, and is given by 

9i = fi0'oi0o,xo,yo)) = -O'o±2p (12) 

xi = x'o (13) 

yi = y'o (14) 

where 6q is the incidence angle, {xi, yi) is the collision point, and (3 is the internal angle of the rombus. The 
± sign depends on the side on which the bounce occours, hence C is piecewise linear in 9'^. 

In order to compute the Lyapunov exponents, we need to evaluate the jacobian matrix J ~ JcJf, product 
of a free-flight Jp and of a collision part Jc, where 

/ de^ 09^ 09^ 
I dOo dxo 8yo 

£ §S I ' Jc = I 1 I . (15) 



Jf 



dy'a dy'o 
36 Q dxQ 




Taking the product of matrices like J at each collision, along an entire trajectory, the Lyapunov exponents 
•^2, A3, A4 of eq.(H) can be computed. The remaining one (corresponding to the conserved kinetic energy) is 
zero. 

Now, consider the relation between the phase space contraction rate and the particle momentum |^,|l6| : 

div(p, q) = -a = -ep^ , (16) 

which holds at all times. This relation implies that the average of div(p, q) (whose asymptotic limit is the 
sum of the Lyapunov exponents) has sign opposite to that of the current and, for both the finite time and 
the asymptotic Lyapunonv exponents, we obtain: 
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Aa;(n) ^ \i{n) + \2{n) + Xajn) + Xjjn) ^^^^ 

Here Ax{n) is the distance travelled in real space in the time t„ corresponding to n collisions, and Xi{n) is 
the i-th finite time exponent. This result is exact, and does not require any conditon on the value of the field. 
It follows that, because a positive current is observed in the stationary state as well as in the transient, the 
sum of the Lyapunonv exponents is always negative. Therefore, one exponent at least is negative. Moreover, 
in the infinite time limit, we know that two exponents vanish, leaving some uncertainty only on the value 
of the remaining exponent. Numerically, we found that the asymptotic value of this Lyapunov exponent is 
either negative or very close to zero. 
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FIG. 1. The modified Ehrenfest gas. Along the paper, the side of the elementary cell is set to 1.291, while the 
semiaxis of the rhombus are chosen to be 1.1 and 0.7573 respectively. 
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FIG. 2. The first 5000 iterates of the bounce map B for e = 0.01 (left panel) starting from a random initial 
condition, and the last 5000 out of a trajectory of 10 collisions (right panel). 
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FIG. 3. The variance cr^ for e = 0.075 and ensemble size A'^ = 300 shows a crossover from diffusive to balhstic 
behaviour at n « ric 
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FIG. 4. The crossover time ric as a function of e. The dashed line corresponds to ric = 2.15 e 



11 



40 



30 



Q. 

V 



1 20 
£ 10 



OT=253.0 160 coll. 
Ot=392.2 240 coll. 
At=522.8 320 coll. 



A_ OA' 



30 



20 - 



10 - 



ox=125.5 80 coll. 
Ox=261 .6 160 coll. 
At=392.4 240 coll. 
+ 1=523.2 320 coll. 





FIG. 5. 

prediction 
field. 



20 

z =<p^>y<p^> 

Test of the fluctuation relation for e 



40 



10 



20 



30 



0.005 (left) and e 



z =<p^>/<p^> 

0.01 (right): the solid line is the theoretical 



Large z-values are more affected by statistical errors but its range of validity of it is larger for smaller 




12 



